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This is a pedagogical review of the recent observational data obtained from 
type la supernova surveys that support the accelerating expansion of the 
universe. The methods for the analysis of the data are reviewed and the the- 
oretical implications obtained from their analysis are discussed. 



1 Introduction 

Recent distance-redshift surveys P El EJ El EJ of cosmologically distant 
Type la supernovae (Snla) have indicated that the universe has recently (at 
\Q , redshift z ~ 0.5) entered a phase of accelerating expansion. This expansion 

has been attributed to a dark energy |7j component with negative pressure 
^ <— | which can induce repulsive gravity and thus cause accelerated expansion. The 

evidence for dark energy has been indirectly verified by Cosmic Microwave 
Background (CMB) 8 and large scale structure |U] observations. 
' The simplest and most obvious candidate for this dark energy is the 

\ cosmological constant ^U] with equation of state w = 2 = — 1. The ex- 

tremely fine tuned value of the cosmological constant required to induce the 
observed accelerated expansion has led to a variety of alternative models 
where the dark energy component varies with time. Many of these mod- 
els make use of a homogeneous, time dependent minimally coupled scalar 
field <j) f quintessence |TT1 H%] ) whose dynamics is determined by a specially de- 
signed potential V{(f>) inducing the appropriate time dependence of the field 
equation of state w(z) — jj^j- Given the observed w(z), the quintessence 
potential can in principle be determined. Other physically motivated mod- 
els predicting late accelerated expansion include modified gravitv fDfl 1141 ITS*] . 
Chaplygin gas^Hj, Cardassian cosmology [T7j. theories with compactified ex- 
tra dimensions |18l IT!?j . braneworld models 2Qj etc. Such cosmological models 
predict specific forms of the Hubble parameter H(z) as a function of redshift 
z. The observational determination of the recent expansion history H(z) is 
therefore important for the identification of the viable cosmological models. 
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The most direct and reliable method to observationally determine the 
recent expansion history of the universe H(z) is to measure the redshift z and 
the apparent luminosity of cosmological distant indicators (standard candles) 
whose absolute luminosity is known. The luminosity distance vs. redshift is 
thus obtained which in turn leads to the Hubble expansion history H{z). 

The goal of this review is to present the methods used to construct the 
recent expansion history H(z) from Snla data and discuss the most recent 
observational results and their theoretical implications. In the next section I 
review the method used to determine H(z) from cosmological distance indi- 
cators and discuss Snla as the most suitable cosmological standard candles. 
In section 3 I show the most recent observational results for H(z) and discuss 
their possible interpretations other than accelerating expansion. In section 4 
I discuss some of the main theoretical implications of the observed H (z) with 
emphasis on the various parametrizations of dark energy (the simplest being 
the cosmological constant). The best fit parametrizations are shown and their 
common features are pointed out. The physical origin of models predicting 
the best fit form of H{z) is discussed in section 5 where I distinguish between 
minimally coupled scalar fields (quintessence) and modified gravity theories. 
An equation of state of dark energy with w < — 1 is obtained by a specific type 
of dark energy called phantom energy |21j . This type of dark energy is faced 
with theoretical challenges related to the stability of the theories that predict 
it. Since however the Snla data are consistent with phantom energy it is in- 
teresting to investigate the implications of such an energy. These implications 
are reviewed in section 6 with emphasis to the Big Rip future singularity im- 
plied by such models as the potential death of the universe. Finally, in section 
7 I review the future observational and theoretical prospects related to the 
investigation of the physical origin of dark energy and summarize the main 
conclusions of this review. 



2 Expansion History from the Luminosity Distances of 
Snla 

Consider a luminous cosmological object emitting at total power L (absolute 
luminosity) in radiation within a particular wavelength band. Consider also 
an observer (see Fig. 1) at a distance dr, from the luminous object. In a static 
cosmological setup, the power radiated by the luminous object is distributed 
in the spherical surface with radius dh and therefore the intensity I (apparent 
luminosity) detected by the observer is 

"A m 

The quantity 
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Fig. 1. The luminosity distance obtaned from the apparent and absolute luminosi- 
ties 



(2) 

is known as the luminosity distance to the luminous object and in a static uni- 
verse it coincides with the actual distance. In an expanding universe however, 
the energy of the radiation detected by the observer has been reduced not 
only because of the distribution of photons on the spherical surface but also 
because the energy of the photons has been redshiftcd while their detection 
rate is reduced compared to their emission rate due to the cosmological ex- 
pansion [22] . Both of these expansion effects give a reduction of the detected 
energy by a factor = (1 + z) where a(t) is the scale factor of the universe 
at cosmic time t and to is the present time. Usually a is normalized so that 
a(io) = 1- Thus the detected apparent luminosity in an expanding background 
may be written as 

1 = 4ira(t ) 2 x{z) 2 (l + z) 2 (3) 

where x(z) is the comoving distance to the luminus object emitting with 
redshift z. This implies that in an expanding universe the luminosity distance 
cIl(z) is related to the comoving distance x(z) by the relation 




<1l{z) = x(z)(l + z) 



(4) 
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Using eq. Q and the fact that light geodesies in a flat expanding background 
obey 

c dt — a(z) dx(z) (5) 

it is straightforward to eliminate x(z) and express the expansion rate of the 
universe H(z) = ~(z) at a redshift z (scale factor a = j-pi) m terms of the 
observable luminosity distance as 

dz 1 + z 

This is an important relation that connects the theoretically predictable Hub- 
ble expansion history H(z) with the observable luminosity distance in 
the context of a spatially flat universe. Therefore, if the absolute luminosity 
of cosmologically distant objects is known and their apparent luminosity is 
measured as a function of redshift, eq. J2J) can be used to calculate their lu- 
minosity distance cZl(z) as a function of redshift. The expansion history H(z) 
can then be deduced by differentiation with respect to the redshift using eq. 
10. Reversely, if a theoretically predicted H(z) is given, the corresponding 
predicted di,{z) is obtained from © by integrating H{z) as 

x f* dz' 

This predicted dh(z) can be compared with the observed dh{z) to test the con- 
sistency of the theoretical model with observations. In practice astronomers 
do not refer to the ratio of absolute over apparent luminosity. Instead they 
use the difference between apparent magnitude m and absolute magnitude M 
which is connected to the above ratio by the relation 

m-M = 2.5 log 10 (j) (8) 

A particularly useful diagram which illustrates the expansion history of 
the Universe is the Hubble diagram. The x-axis of a Hubble diagram (see 
Fig. 2) shows the redshift z of cosmological luminous objects while the y- 
axis shows the physical distance Ar to these objects. In the context of a 
cosmological setup the redshift z is connected to the scale factor a(t) at the 
time of emission of radiation by 1 + z — ^pj- where to is the present time. On 
the other hand, the distance to the luminous object is related to the time in 
the past t past when the radiation emission was made. Therefore, the Hubble 
diagram contains information about the time dependence of the scale factor 
ait). The slope of this diagram at a given redshift denotes the inverse of the 
expansion rate ~{z) = H(z) ie 
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Fig. 2. The Hubble diagram. In an accelerating universe luminous objects at a given 
rcdshift appear to be dimmer. 



In an accelerating universe the expansion rate H(z) was smaller in the past 
(high rcdshift) and therefore the slope H~ x of the Hubble diagram is larger at 
high rcdshift. Thus, at given redshift, luminous objects appear to be further 
away (dimmer) compared to an empty universe expanding with a constant 
rate (see Fig. 2). 

The luminous objects used in the construction of the Hubble diagram are 
objects whose absolute luminosity is known and therefore their distance can 
be evaluated from their apparent luminosity along the lines discussed above. 
Such objects are known as distance indicators or standard candles. A list of 
common distance indicators used in astrophysics and cosmology is shown in 
Table 1 along with the range of distances where these objects are visible and 
the corresponding accuracy in the determination of their absolute magnitude. 
As shown in Table 1 the best choice distance indicators for cosmology are 
Snla not only because they are extremely luminous (at their peak they are as 
luminous as a bright galaxy) but also because their absolute magnitude can 
be determined at a high accuracy. 

Type la supernovae emerge in binary star systems where one of the com- 
panion stars has a mass below the Chandrasekhar limit IAMq and therefore 
ends up (after hydrogen and helium burning) as white dwarf supported by 
degeneracy pressure. Once the other companion reaches its red giant phase 
the white dwarf begins gravitational striping of the outer envelop of the red 
giant thus accreting matter from the companion star. Once the white dwarf 
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Table 1. Extragalactic distance indicators (from Ref. |23|) 



Technique 


Range of distance 


Accuracy (la) 


Cepheids 


< LMC to 25 Mpc 


0.15 mag 


SNIa 


4 Mpc to > 2 Gpc 


0.2 


mag 


Expand. Phot. Meth./Snll 


LMC to 200 Mpc 


0.4 


mag 


Planetary Nebulae 


LMC to 20 Mpc 


0.1 


mag 


Surf. Brightness Fluct 


1 Mpc to 100 Mpc 


0.1 


mag 


Tully Fisher 


1 Mpc to 100 Mpc 


0.3 


mag 


Brightest Cluster Gal. 


50 Mpc to 1 Gpc 


0.3 


mag 


Glob. Cluster Lum. Fun. 


1 Mpc to 100 Mpc 


0.4 


mag 


Sunyaev-Zeldovich 


100 Mpc to > 1 Gpc 


0.4 


mag 


Gravitational Lensing 


5 Gpc 


0.4 


mag 



reaches a mass equal to the Chandrasekhar limit, the degeneracy pressure is 
unable to support the gravitational pressure, the white dwarf shrinks and in- 
creases its temperature igniting carbon fussion. This leads to violent explosion 
which is detected by a light curve which rapidly increases luminosity in a time 
scale of less than a month, reaches a maximum and disappears in a timescale 
of 1-2 months (see Fig. 3). Type la are the preferred distance indicators for 



-20 r 




Fig. 3. Typical Snia light-curve. 



cosmology for several reasons: 
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1. They are exceedingly luminous. At their peak luminosity they reach an 
absolute magnitude of M ~ — 19 which corresponds to about 10 10 Af Q . 

2. They have a relatively small dispersion of peak absolute magnitude. 

3. Their explosion mechanism is fairly uniform and well understood. 

4. There is no cosmic evolution of their explosion mechanism according to 
known physics. 

5. There are several local Snla to be used for testing Snla physics and for 
calibrating the absolute magnitude of distant Snla. 

On the other hand, the main problem for using Snla as standard candles 
is that they are not easy to detect and it is impossible to predict a Snla 
explosion. In fact the expected number of Snla exploding per galaxy is 1-2 
per millenium. It is therefore important to develop a search strategy in order 
to efficiently search for Snla at an early stage of their light curve. The method 
used (with minor variations) to discover and follow up photometrically and 
spectroscopically Snla consists of the following steps ^ |3 |3J : 

1 . Observe a number of wide fields of apparently empty sky out of the plane 
of our Galaxy. Tens of thousands of galaxies are observed in a few patches 
of sky. 

2. Come back three weeks later (next new moon) to observe the same galaxies 
over again. 

3. Subtract images to identify on average 12-14 Snla. 

4. Schedule in advance follow up photometry and spectroscopy on these Snla 
as they brighten to peak and fade away. 

Given the relatively short time difference (three weeks) between first and sec- 
ond observation, most Snla do not have time to reach peak brightness so al- 
most all the discoveries are pre-maximum. This strategy turns a rare, random 
event into something that can be studied in a systematic way. This strategy 
is illustrated in Figs 4 and 5 (from Ref. [21] )• The outcome of this observation 
strategy is a set of Snla light curves in various bands of the spectrum (see 
Fig. 6). These light curves are very similar to each other and their peak ap- 
parent luminosity could be used to construct the Hubble diagram assuming a 
common absolute luminosity. 

Before this is done however a few corrections must be made to take into 
account the minor intrinsic absolute luminosity differences (due to composi- 
tion differences) among Snla as well as the radiation extinction due to the 
intergalactic medium. Using samples of closeby Snla it has been empirically 
observed that the minor differences of Snla absolute luminosity are connected 
with differences in the shape of their light curves. Broad slowly declining 
light curves (stretch factor s > 1) correspond to brighter Snla while narrower 
rapidly declining light curves (stretch factor s < 1) correspond to intrinsically 
fainter Snla. This stretch factor dependence of the Snla absolute luminosity 
has been verified using closeby Snla It was shown that contraction of 
broad light curves while reducing peak luminosity and stretching narrow light 
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Fig. 5. Supernova 1998ba, an example of a supernova discovery using the search 
strategy described in the text involving subtraction of images. 
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Fig. 6. A set of light curves from SN2001el in various bands of the spectrum. 

curves while increasing peak luminosity makes these light curves coincide (see 
Fig. 7). In addition to the stretch factor correction an additional correction 
must be made in order to compare the light curves of high redshift Snla with 
those of lower redshift. In particular all light curves must be transformed to 
the same reference frame and in particular the rest frame of the Snla. For 
example a low redshift light curve of the blue B band of the spectrum should 
be compared with the appropriate red R band light curve of a high redshift 
Snla. The transformation also includes correction for the cosmic time dilation 
(events at redshift z last 1 + z times longer than events at z ~ 0). These cor- 
rections consist the K- correction and is used in addition to the stretch factor 
correction discussed above. The K-correction transformation is illustrated in 
Fig.8. 



3 Observational Results 

The first project in which Snla were used to determine the cosmological 
constant energy was the research from Perlmutter et al. in 1997 26 . The 
project was known as the Supernova Cosmology Project (SCP). Applying 
the above described methods they discovered seven distant Snla at redshift 
0.35 < z < 0.65. When discovered, the supernovae were followed for a year 
by different telescopes on earth to obtain good photometry data in differ- 
ent bands, in order to measure good magnitudes. The Hubble diagram they 
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Fig. 7. Left:The range of lightcurve for low-redshift supernovae discovered by the 
Calan/Tololo Supernova Survey. At these redshifts, the relative distances can be 
determined (from redshift), so their relative brightnesses are known. Right: The 
same lightcurves after calibrating the supernova brightness using the stretch of the 
timescale of the lightcurve as an indicator of brightness (and the color at peak as 
an indicator of dust absorption) 



constructed was consistent with standard Friedman cosmology without dark 
energy or cosmological constant. 

A year after their first publication, Perlmutter et al. published in Nature 
an update on their initial results. They had included the measurements of a 
very high-redshifted z — 0.83 Supernova la. This dramatically changed their 
conlusions. The standard decelerating Friedman cosmology was rulled out 
at about 99% confidence level. The newly discovered Supernova indicated a 
universe with accelerating expansion dominated by dark energy. These results 
were confirmed independently by another pioneer group (High-z Supernova 
Search Team (HSST)) searching for Snla and measuring the expansion history 
H(z) (Riess et al. in 1998 0). They had discovered 16 Snla at 0.16 < z < 0.62 
and their H(z) also indicated accelerating expansion ruling out for a flat 
universe. Their data also permitted them to definitely rule out decelerating 
Friedman cosmology at about 99% confidence level. 

In 2003 Tonry et al. [3] reported the results of their observations of eight 
newly discovered Snla. These Snla were found in the region 0.3 < z < 1.2. 
Together with previously acquired Snla data they had a data set of more than 
100 Snla. This dataset confirmed the previous findings of accelerated expan- 
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Fig. 8. Slightly different parts of the supernova spectrum are observed through 
the B filter transmission function at low redshift (upper panel) and through the R 
filter transmission function at high redshift (lower panel). This small difference is 
accounted for by the 'cross-filter K-correction' 25: . 



sion and gave the first hints of decelerated expansion at redshifts z ~ 0.6 
when matter is expected to begin dominating over dark energy. This transi- 
tion from decelerating to accelerating expansion was confirmed and pinpointed 
accurately by Riess et al. in 2004 5 who included in the analysis 16 new high- 
redshift Snla obtained with HST and reanalyzed all the available data in a 
uniform and robust manner constructing a robust and reliable dataset consist- 
ing of 157 points known as the Gold dataset. These Snla included 6 of the 7 
highest redshift Snla known with z > 1.25. With these new observations, they 
could clearly identify the transition from a decelerating towards an accelerat- 
ing universe to be at z — 0.46 ±0.13. It was also possible to rule out the effect 
of dust on the dimming of distant Snla, since the accelerating/decelerating 
transition makes the effect of dimming inverse. The Hubble diagram obtained 
from the Gold dataset is shown in Fig. 9 where the corrected apparent mag- 
nitude m(z) of the 157 Snla is plotted versus the redshift z. The apparent 
magnitude m(z) is related to the corresponding luminosity distance dh of the 
Snla by 

m(z) =M + 5Zo0 lo [%^] + 25 (10) 
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Fig. 9. The apparent magnitude m(z) vs redshift as obtained from the Gold dataset. 
It is not easy to distinguish between accelerating and decelerating expansion in such 
a diagram. 



where M is the absolute magnitude which is assumed to be constant for 
standard candles like Snla after the corrections discussed in section 2 are 
implemented. 

A potential problem of plots like the one of Fig. 9 is that it is not easy 
to tell immediately if the data favor an accelerating or decelerating universe. 
This would be easy to tell in the Hubble diagram of Fig. 2 where the distance 
is plotted vs redshift and is superposed with the distance-redshift relation 
dT npv (z) of an empty universe with H{z) constant. An even more efficient 
plot for such a purpose would be the plot of the ratio J^}u), (or its log w ) 
which can immediately distinguish accelerating from decelerating expansion 
by comparing with the jt^h), = 1 line. Such a plot is shown in Fig. 10 |S] 

using both the raw Gold sample data and the same data binned in redshift 
bins. 

The lines of zero acceleration, constant acceleration and constant decel- 
eration are also shown for comparison. Clearly the best fit is obtained by an 
expansion which is accelerating at recent times (z ~ 0.5) and decelerating at 

earlier times (z ~ 0.5) when matter is expected to dominate. 

The interpretation of the data assuming that the observed dimming at high 
redshift is due to larger distance may not be the only possible interpretation. 
The most natural alternative interpretations however have been shown to lead 
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Fig. 10. The reduced Hubble diagram used to distinguish between accelerating and 
decelerating expansion 5 . 



to inconsistencies and none of them has been favored as a viable alternative 
at present. These alternative interpretations include the following: 

• Intergalactic Dust: Ordinary astrophysical dust does not obscure equally 
at all wavelengths, but scatters blue light preferentially, leading to the well- 
known phenomenon of "reddening" . Spectral measurements |S] reveal a 
negligible amount of reddening, implying that any hypothetical dust must 
be a novel "grey" variety inducing no spectral distortions |27| . 

• Grey Dust: Grey dust is highly constrained by observations: first, it pre- 
dicts further increase of dimming at higher redshifts z ~ 0.5 which is not 
observed; and second, intergalactic dust would absorb ultraviolet/optical 
radiation and re-emit it at far infrared wavelengths, leading to stringent 
constraints from observations of the cosmological far-infrared background. 
Thus, while the possibility of obscuration has not been entirely eliminated, 
it requires a novel kind of dust which is already highly constrained (and 
may be convincingly ruled out by further observations). 

• Evolution of Snla: The supernova search teams have found consistency 
in the spectral and photometric properties of Snla over a variety of red- 
shifts and environments 5 (e.g., in elliptical vs. spiral galaxies). Thus 
despite the relevant tests there is curently no evidence that the observed 
dimming can be attributed to evolution of Snla. 
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According to the best of our current understanding, the supernova results 
indicating an accelerating universe seem likely to be trustworthy. Needless 
to say, however, the possibility of a neglected systematic effect can not be 
definitively excluded. Future experiments, discussed in section 7 will both 
help us improve our understanding of the physics of supernovae and allow a 
determination of the distance/redshift relation to sufficient precision to dis- 
tinguish between the effects of an accelerating universe and those of possible 
astrophysical phenomena. 

4 Dark Energy and Negative Pressure 

Our current knowledge of the expansion history of the universe can be summa- 
rized as follows: The universe originated at an initial state that was very close 
to a density singularity known as the Big Bang. Soon after that it entered 
a phase of superluminal accelerating expansion known as inflation. During 
inflation causally connected regions of the universe exited out of the horizon, 
the universe approached spatial flatness and the primordial fluctuations that 
gave rise to structure were generated. At the end of inflation the universe 
was initially dominated by radiation and later by matter whose attractive 
gravitational properties induced a decelerating expansion. 

The Snla data discussed in section 3 (along with other less direct cos- 
mological observations HHU) strongly suggest that the universe has recently 
entered a phase of accelerating expansion at a redshift z ~ 0.5. This accelerat- 
ing expansion can not be supported by the attractive gravitational properties 
of regular matter. The obvious question to address is therefore 'What are 
the properties of the additional component required to support this acceler- 
ation?'. To address this question we must consider the dynamical equation 
that determines the evolution of the scale factor a(t). This equation is the 
Friedman equation which is obtained by combining General Relativity with 
the cosmological principle of homogeneity and isotropy of the universe, ft may 
be written as 



where pi and pi are the densities and pressures of the contents of the universe 
assumed to behave as ideal fluids. The only directly detected fluids in the uni- 
verse are matter (p m ,Pm = 0) and the subdominant radiation (p r ,Pr = Pr/3). 
Both of these fluids are unable to cancel the minus sign on the rhs of the 
Friedman equation and can therefore only lead to decelerating expansion. Ac- 
celerating expansion in the context of general relativity can only be obtained 
by assuming the existence of an additional component {px,Px — wpx) termed 
'dark energy' which could potentially change the minus sign of eq. (|11|1 and 
thus lead to accelerating expansion. Assuming a positive energy density for 




(11) 
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dark energy (required to achieve flatness) it becomes clear that negative pres- 
sure is required for accelerating expansion. In fact, writing the Friedman eq. 
I|ll|) in terms of the dark energy equation of state parameter w as 

a iirG . , . , , 

- = —[ Pm + Px{l + M\ (12) 

a 3 

it becomes clear that a w < — i is required for accelerating expansion implying 
repulsive gravitational properties for dark energy. 

The redshift dependence of the dark energy can be easily connected to 
the equation of state parameter w by combining the energy conservation 
d(pxa 3 ) — —p x d(a 3 ) with the equation of state px — wpx as 

px ~ a- 3 ^ = (1 + z) 3 < 1+ ™) (13) 

This redshift dependence is related to the observable expansion history H(z) 
through the Friedman equation 

2 = = 8*G oo )3 + px[a)] = H 2 [[2om{l + z)3 + Qx{z)] (14) 

where the density parameter Q = p ? - for matter is constrained by large 
scale structure observations to a value (prior) i7om — 0.3. Using this prior, 
the dark energy density parameter f2x{z) = and the corresponding 

equation of state parameter w may be constrained from the observed H(z). 

In addition to Slx(z), the luminosity distance- redshift relation (Il(z) ob- 
tained from Snla observations can constrain other cosmological parameters. 
The only parameter however obtained directly from cIl(z) (using eq. ©) is 
the Hubble parameter H(z). Other cosmological parameters can be obtained 
from H(z) as follows: 

• The age of the universe io is obtained as: 

f°° dz 
Jo { l + z )H{ z ) 

• The present Hubble parameter Ho = H(z = 0). 

• The deceleration parameter q(z) = 

q(z = {1 + z)— 1 (16) 

dz 

and its present value go = l( z = 0). 

• The density parameters for matter and dark energy are related to H[z) 
through the Friedman equation Q14JI. 

• The equation of state parameter w(z) obtained as |281 129| 

W (z) - P -^l = + (17) 

1 ' p X (z) l_(Sl)2^ 0m (l + z)3 V > 

obtained using the Friedman equations l|12f> and (|14fl . 
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The most interesting parameter from the theoretical point of view (apart 
from H(z) itself) is the dark energy equation of state parameter w(z). This 
parameter probes directly the gravitational properties of dark energy which 
are predicted by theoretical models. The downside of it is that it requires two 
differentiations of the observable <1l(z) to be obtained and is therefore very 
sensitive to observational errors. 

The simplest form of dark energy corresponds to a time independent en- 
ergy density obtained when w = — 1 (see eq. (f 13|) ^ - This is the well known 
cosmological constant which was first introduced by Einstein in 1917 two 
years after the publication of the General Relativity (GR) equation 

= kT^ (18) 

where k = 8ttG/c 2 . At the time the 'standard' cosmological model was a 
static universe because the observed stars of the Milky Way were found to 
have negligible velocities. The goal of Einstein was to apply GR in cosmology 
and obtain a static universe using matter only. It became clear that the attrac- 
tive gravitational properties of matter made it impossible to obtain a static 
cosmology from (|18f) . A repulsive component was required and at the time of 
major revolutions in the forms of physical laws it seemed more natural to ob- 
tain it by modifying the gravitational law than by adding new forms of energy 
density. The simplest generalization of eq. (|18fl involves the introduction of a 
term proportional to the metric g M „. The GR equation becomes 

- Ag^ = kT^ v (19) 

where A is the cosmological constant. The repulsive nature of the cosmological 
constant becomes clear by the metric of a point mass (Schwarschild-de Sitter 
metric) which, in the Newtonian limit leads to a gravitational potential 

(20, 

r b 

which in addition to the usual attractive gravitational term has a repulsive 
term proportional to the cosmological constant A. This repulsive gravitational 
force can lead to a static (but unstable) universe in a cosmological setup 
and in the presence of a matter fluid. A few years after the introduction 
of the cosmological constant by Einstein came Hubble's discovery that the 
universe is expanding and it became clear that the cosmological constant was 
an unnecessary complication of GR. It was then that Einstein (according to 
Gamow's autobiography) called the introduction of the cosmological constant 
'the biggest blunder of my life'. In a letter to Lemaitre in 1947 Einstein wrote: 
'Since I introduced this term I had always had a bad conscience. I am unable 
to believe that such an ugly thing is actually realized in nature'. As discussed 
below, there is better reason than ever before to believe that the cosmological 
constant may be non-zero, and Einstein may not have blundered after all. 
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If the cosmological constant is moved to the right hand side of eq. I|19fl it 
may be incorporated in the energy momentum tensor as an ideal fluid with 
Pa = and w = — 1. In the context of field theory such an energy momen- 
tum tensor is obtained by a scalar field <fi with potential V(4>) at its vacuum 
state 4>o ie d^cj) = and = — y(0 o )g M „. Even though the cosmologi- 
cal constant may be physically motivated in the context of field theory and 
consistent with cosmological observation there are two important problems 
associated with it: 

• Why is it so incredibly small? Observationally, the cosmological constant 
density is 120 orders of magnitude smaller than the energy density associ- 
ated with the Planck scale - the obvious cut off. Furthermore, the standard 
model of cosmology posits that very early on the universe experienced a 
period of inflation: A brief period of very rapid acceleration, during which 
the Hubble constant was about 52 orders of magnitude larger than the 
value observed today How could the cosmological constant have been so 
large then, and so small now? This is sometimes called the cosmological 
constant problem. 

• The 'coincidence problem': Why is the energy density of matter nearly 
equal to the dark energy density today? 

Despite the above problems and given that the cosmological constant is the 
simplest dark energy model, it is important to investigate the degree to which 
it is consistent with the Snla data. I will now describe the main steps involved 
in this analysis. According to the Friedman equation the predicted Hubble 
expansion in a flat universe and in the presence of matter and a cosmological 
constant is 

= oj = 8^7 ao )3 + A H 2 [Qom{l + z) 3 + Qa] (21) 
a o a 6 

where Qa = PA and 

tfom + A = 1 (22) 

This is the LCDM (yl+Cold Dark Matter) which is currently the minimal 
standard model of cosmology. The predicted H(z) has a single free parameter 
which we wish to constrain by fitting to the Snla luminosity distance-redshift 
data. 

Observations measure the apparent luminosity vs redshift l(z) or equiv- 
alently the apparent magnitude vs redshift m(z) which are related to the 
luminosity distance by 

2.5Zo 3l0 ( A) = m{z) — M — 25 = 5log 10 (^^) (23) 
lyz) Mpc 

From the theory point of view the predicted observable is the Hubble param- 
eter (|21ll which is related to the theoretically predicted luminosity distance 
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dh( z ) by eq. (7J. In this case d^z) depends on the single parameter f2 0m and 
takes the form 

f z dz' 

d L (z;n 0m ) th =c(l+z) - (24) 

JO " \ Z I l2 0rn) 

Constraints on the parameter J?om are obtained by the maximum likelihood 
method 30 which involves the minimization of the x 2 (^0m) defined as 

2tn x [dL(z)obs - d L (z; f2 0m ) th } 2 

X (MOm) — 2 \ Zb > 

i=l « 

where N is the number of the observed Snla luminosity distances and Oi are 
the corresponding la errors which include errors due to flux uncertainties, in- 
ternal dispersion of Snla absolute magnitude and peculiar velocity dispersion. 
If flatness is not imposed as a prior through eq. I|22|l then dL{z)th depends 
on two parameters (f2o m and £2 A ) and the relation between dh(z\ fiom, & A )th 
and H(z] J?omj &a) takes the form 

d L (z)th = -— C(1 + ^ = sin[y/f2 0m + TT A - 1 / dz'-^-A (26) 
V^om + n A -i Jo H ( z ) 

In this case the minimization of eq. (|25|l leads to constraints on both f2o m and 
Qa- This is the only direct and precise observational probe that can place 
constraints directly on Q A - Most other observational probes based on large 
scale structure observations place constraints on f2o m which are indirectly 
related to fl A in the context of a flatness prior. 

As discussed in section 2 the acceleration of the universe has been con- 
firmed using the above maximum likelihood method since 1998 ^IH- Even 
the early datasets of 1998 ^ |2] were able to rule out the flat matter domi- 
nated universe (SCDM: f] 0m = 1, fl A = 0) at 99% confidence level. The latest 
datasets are the Gold dataset (N = 157 in the redshift range < z < 1.75) 
discussed in section 2 and the first year SNLS (Supernova Legacy Survey) 
dataset which consists of 71 datapoints in the range < z < 1 plus 44 pre- 
viously published closeby Snla. The 68% and 95% \ 2 contours in the (Qa m 
and £2 a) parameter space obtained using the maximum likelihood method 
are shown in Fig. 11 for the SNLS dataset, a truncated version of the Gold 
dataset (TG) with < z < 1 and the Full Gold (FG) dataset. The following 
comments can be made on these plots: 

• The two versions of the Gold dataset favor a closed universe instead of a 
flat universe {Qfg = 2.16 ± 0.59, f2f a f = 1.44 ± 0.44). This trend is not 
realized by the SNLS dataset which gives flf£ LS = !-07 ± 0.52. 

• The point corresponding to SCDM (J?om, ft a) = (L 0) is ruled out by all 
datasets at a confidence level more than 10cr. 

• If we use a prior constraint of flatness fio m + Q A = 1 thus restricting on 
the corresponding dotted line of Fig. 1 and using the parametrization 
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Fig. 11. The 68% and 95% \ 2 contours in the (Oo m and Qa) parameter space 
obtained using the SNLS, TG and FG datasets (from Ref. |5T)L 



H(z) 2 = i/ 2 [r2 0m (l + 2) 2 + (1 - flom)] (27) 
we find minimizing x 2 (^0m) of eq i|25[) 

^om i5 = °- 26 ± °- 04 ( 28 ) 
/2™ = 0.30 ± 0.05 (29) 
f2£G = 0.31 ± 0.04 (30) 

These values of J?om are consistent with corresponding constraints from 
the CMBJ5] and large scale structure observations [D]. 

Even though LCDM is the simplest dark energy model and is currently con- 
sistent with all cosmological observations (especially with the SNLS dataset) 
the question that may still be address is the following: 'Is it possible to get 
better fits (lowering \ 2 further) with different H(z) parametrizations and if 
yes what are the common features of there better fits?' The strategy towards 
addressing this question involves the following steps: 

• Consider a physical model and extract the predicted recent expansion his- 
tory H(z; ax, a,2, a n ) as a function of the model parameters ai, a-2, ■ «n- 
Alternatively a model independent parametrization for H{z; ax, o, n ) 
(or equivalently w(z; ax, fl2, a n )) may be constructed aiming at the best 
possible fit to the data with a small number of parameters (usually 3 or 
less). 

• Use eq. @ to obtain the theoretically predicted luminosity distance 
dh(z\ ax, Q,2, a n )th as a function of z. 

• Use the observed luminosity distances dh{zi)obs to construct \ 2 along 
the lines of eq. I|25|) and minimize it with respect to the parameters 
ax, a.2, a n . 
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• From the resulting best fit parameter values ai,a 2 ,...,a„ (and their er- 
ror bars) construct the best fit H (z; Oi, 0,2, a n ), <1l(z; a 1; a 2 , a„) and 
u>(z; 5i, 02, a n ). The quality of fit is measured by the depth of the min- 
imum of x 2 ie xL„(ai,a 2 , -,a n )- 

Most useful parametrizations reduce to LCDM of eq. I|21[l for specific param- 
eter values giving a Xlcdm f° r these parameter values. Let 

Axlcdm = Xmm(«ii 5 2, — , a n ) - Xlcdm (31) 

The value of Ax\ CDM is usually negative since % 2 is usually further reduced 
due to the larger number of parameters compared to LCDM. For a given 
number of parameters the value of Ax\ CDM gives a measure of the probability 
of having LCDM physically realized in the context of a given parametrization 
|32| . The smaller this probability is, the more 'superior' this parametrization 
is compared to LCDM. For example for a two parameter parametrization and 
\Ax 2 LCDM \ > 2.3 the parameters of LCDM are more than 1<t away from the 
best fit parameter values of the given parametrization. This statistical test 
has been quantified in Ref. |32| and applied to several H (z) parametrizations. 

As an example let us consider the two parameter polynomial parametriza- 
tion allowing for dark energy evolution 

H(zf = Hi [f2 0m (l + zf + a 2 (l + zf + o x (l + z) + (1 - a 2 - a x - f2 0m )} (32) 

in the context of the Full Gold dataset. Applying the above described x 2 
minimization leads to the best fit parameter values a\ = 1.67 ± 1.03 and 
a 2 = — 4.16± 2.53. The corresponding \Ax\cdm\ 1S found to be 2.9 which 
implies that the LCDM parameters values (a\ — a-i = 0) are in the range of 
1<7 — 2a away from the best fit values. 

The same analysis can be repeated for various different parametrizations in 
an effort to identify the common features of the best fit parametrizations. For 
example two other dynamical dark energy parametrizations used commonly 
in the literature are defined in terms of w(z) as 

• Parametrization A: 

w(z) = wq + wi z (33) 
H 2 (z) = Hllf2 Qm (l + zf + (34) 
+ (1 - Q 0m )(l + z )^ 1+wa - w ^e 3wiz ] 

• Parametrization B: 

w(z) = Wq + Wi— — (35) 
_L i Z 

H 2 {z)=Hl[Q Gm {l + zf + 

+ (1 - n 0m )(l + z )3(l+toa+«i) e 3iw[l/(l+*)-l]] ( 36 ) 
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Fig. 12. The best fit forms of w(z) obtained from a variety of parametrizations 32 
in the context of the Full Gold dataset. Notice that they all cross the line w = — 1 
also known as the Phantom Divide Line (PDL). 

where the corresponding forms of H(z) are derived using eq. <|17|1 . The best 
fit forms of w{z) obtained from a variety of these and other parametrizations 
|3*2*| in the context of the Full Gold dataset are shown in Fig. 12. Even though 
these best fit forms appear very different at redshifts z > 0.5 (mainly due 
to the two derivatives involved in obtaining w(z) from cLl{z)), in the range 
< z < 0.5 they appear to have an interesting common feature: they all cross 
the line w = — 1 also known as the Phantom Divide Line (PDL). As discussed 
in the next section this feature is difficult to reproduce in most theoretical 
models based on minimally coupled scalar fields and therefore if it persisted 
in other independent datasets it could be a very useful tool in discriminating 
among theoretical models. Unfortunately if the same analysis is made in the 
context of the more recent SNLS dataset it seems that this common feature 
does not persist. In Fig. 13 the best fit w(z) (along with the la error region) 
is shown in the context of three different datasets (in analogy with Fig. 11) for 
the there different parametrizations (A, B and polynomial of eq. i|32|) (called 
C in Fig. 13)). Even though the crossing of the PDL is realized at best fit for 
both the FG and TG datasets it is not realized at best fit when the SNLS is 
used. Thus we must wait until further Snla datasets are released before the 
issue is settled. In Fig. 14 I show the la and 2a \ 2 contours corresponding 
to parametrizations A and B with a prior of J?om = 0.24 confirming the fact 
that the SNLS dataset provides best fit parameter values that are almost 
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Fig. 13. The best fit w(z) (along with the lcr error (shaded region)) is shown in 
the context of three different datasets (in analogy with Fig. 11) for there different 
parametrizations (A, B and C) |31) . 

identical to those corresponding to LCDM (wq = —1, w\ = 0) despite the 
dynamical degrees of freedom incorporated in the parametrizations A and B. 
It should be pointed out however that despite the differences in the best fit 
parametrizations, the three datasets (SNLS, TG and FG) are consistent with 
each other at the 95% confidence range (see e.g. Fig. 14) and they are all 
consistent with flat LCDM with J?Q m — 0.3. 

5 Theoretical Models for Dark Energy 

Even though LCDM is the simplest model consistent with current cosmo- 
logical data it is plagued with theoretical fine tuning problems discussed in 
the previous section (the 'coincidence' and the 'cosmological constant' prob- 
lems). In additions dynamical dark energy parametrizations of H(z) provide 
in certain cases significantly better fits to the Snla data. Therefore the inves- 
tigation of physically motivated models that predict a dynamical evolution of 
dark energy is an interesting and challenging problem. 

The role of dark energy can be played by any physical field with posi- 
tive energy and negative pressure which violates the strong energy condition 
p + 3p> (w > — |). Quintessence scalar fields [33] with small positive kinetic 
term (—1 < w < — i) violate the strong energy condition but not the domi- 




Fig. 14. The la and 2a \ 2 contours corresponding to parametrizations A and B 
with a prior of Qom = 0.24. Notice that the SNLS dataset provides best fit parameter 
values that are almost identical to those corresponding to LCDM (wo = — 1, Wi = 0). 



nant energy condition p + p > 0. Their energy density scales down with the 
cosmic expansion and so does the cosmic acceleration rate. Phantom fields 34 
with negative kinetic term (w < —1) violate the strong energy condition, the 
dominant energy condition and maybe physically unstable. However, they 
are also consistent with current cosmological data and according to recent 
studies jHHl EHlliSl they maybe favored over their quintessence counterparts. 

Homogeneous quintessence or phantom scalar fields are described by La- 
grangians of the form 

C = ±^ 2 -V(4>) (37) 

where the upper (lower) sign corresponds to a quintessence (phantom) field 
in equation (|37|l and in what follows. The corresponding equation of state 
parameter is 

P ±\4> 2 -v(^) 

w = - = — 38) 

p ±^ 2 + y(0) 

For quintessence (phantom) models with V(<fi) > (V((f>) < 0) the parameter 
w remains in the range —1 < w < 1. For an arbitrary sign of V((j>) the above 
restriction does not apply but it is still impossible for w to cross the PDL 
it) = —1 in a continous manner. The reason is that for «i = -la zero kinetic 
term ±(f> 2 is required and the continous transition from w < — ltou>>— f (or 
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vice versa) would require a change of sign of the kinetic term. The sign of this 
term however is fixed in both quintessence and phantom models. This difficulty 
in crossing the PDL w = — 1 could play an important role in identifying the 
correct model for dark energy in view of the fact that data favor w ~ — 1 and 
furthermore parametrizations of w(z) where the PDL is crossed appear to be 
favored over the cosmological constant w = — 1 according to the Gold dataset 
as discussed in the previous section. 

It is therefore interesting to consider the available quintessence and phan- 
tom scalar field models and compare the consistency with data of the predicted 
forms of w(z) among themselves and with arbitrary parametrizations of w{z) 
that cross the PDL. This task has been recently undertaken by several authors 
in the context of testing the predictions of phantom and quintessence scalar 

field models 123 ESI- 

As an example we may consider a particular class of scalar field potentials 
of the form 

V{<t>) = s <j> (39) 

where I have followed Ref. EZ] and set <j) — at V = 0. As discussed in 
section 2 (see also Ref. EZj) the field may be assumed to be frozen (0 = 0) 
at early times due to the large cosmic friction H(t). It has been argued '38, 
that such a potential is favored by anthropic principle considerations because 
galaxy formation is possible only in regions where V(<f)) is in a narrow range 
around V = and in such a range any potential is well approximated by a 
linear function. In addition such a potential can provide a potential solution 
to the cosmic coincidence problem [55] . 

The cosmological evolution in the context of such a model 011] is obtained 
by solving the coupled Friedman-Robertson- Walker (FRW) and the scalar 
field equation 

a -r- 1 r& i a\ Q>™ H o f An \ 
a =T 3Ml {(t> +S<t)) -^T- (4°) 

+ 3-0-s = O (41) 
a 

where M p = (8nGy 1/2 is the Planck mass and I have assumed a potential of 
the form 

V{<t>) =Ts<f> (42) 

where the upper (lower) sign corresponds to quintessence (phantom) models. 
The solution of the system (|4U|l - <|4ip for both positive and negative values of 
the single parameter of the model s, is a straightforward numerical problem 
40 which leads to the predicted forms of H(z;s) and w(z;s). These forms 
may then be fit to the Snla datasets for the determination of the best fit 
value of the parameter s. This task has been undertaken in Ref. 001 using 
the Full Gold dataset. The best fit value of s was found to be practically 
indistinguishable from zero which corresponds to the cosmological constant 
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for both the quintessence and the phantom cases. The predicted forms of 
w(z) for a phantom and a quintessence case and s ~ 2 is shown in Fig. 
15. The value of ^Xlcdm ^ s positive in both cases which implies that the 




Fig. 15. The predicted forms of w(z) for a phantom and a quintessence case and 
s ~ 2 provide worse fits to the Gold dataset than LCDM and even worse compared 
to best fit parametrizations that cross the PDL |4U|. 



fit is worse compared to LCDM. The main reason for this is that both the 
quintessence and phantom minimally coupled scalar field models do not allow 
for crossing of the PDL line for any parameter value as discussed above. In 
contrast, the best fit w(z) parametrizations A and B of eqs. (|33[l - i|35|) which 
allow for PDL crossing have a negative A\ 1 LCDM in the context of the Gold 
dataset as shown in Fig. 15 and therefore provide better fits than the field 
theory models. It should be stressed however that in the context of the SNLS 
dataset, parametrizations that allow for crossing of the PDL do not seem to 
have a similar advantage as discussed in the previous section. 

The difficulty in crossing the PDL w = — 1 described above could play 
an important role in identifying the correct model for dark energy in view of 
the fact that data favor w ~ — 1 and furthermore parametrizations of w(z) 
where the PDL is crossed appear to be favored over the cosmological constant 
w = — 1 in the context of the Gold dataset. Even for generalized k-essence 
Lagrangians^J 021 of a minimally coupled scalar field eg 
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£ = \im 2 - V{cj>) (43) 



it has been shown 02j to be impossible to obtain crossing of the PDL. Multiple 
field Lagrangians (combinations of phantom with quintessence fields ^] 23 
0BJ07j) have been shown to in principle achieve PDL crossing but such models 
are complicated and without clear physical motivation (but see |48| for an 
interesting physically motivated model). 

The obvious class of theories that could lead to a solution of the above 
described problem is the non-minimally coupled scalar fields. Such theories 
are realized in a universe where gravity is described by a scalar-tensor theory 
and their study is well motivated for two reasons: 

1. A scalar-tensor theory of gravity is predicted by all fundamental quan- 
tum theories that involve extra dimensions. Such are all known theories 
that attempt to unify gravity with the other interactions (eg supergravity 
(SUGRA), M-theory etc). 

2. Scalar fields emerging from scalar tensor theories (extended quintessence) 
can predict an expansion rate H (z) that violates the inequality 

az 

which is equivalent to crossing the PDL w = — 1 (see eg Ref. |49j'). 

In fact it has been shown in Ref. |49| that in contrast to minimally coupled 
quintessence, scalar tensor theories can reproduce the main features of the 
best fit Hubble expansion history obtained from the Gold dataset. However, 
the precise determination of the scalar tensor theory potentials requires more 
accurate Snla data and additional cosmological observational input. 



6 The Fate of a Phantom Dominated Universe: Big Rip 

As discussed in section 4 the Gold dataset favors a dynamical dark energy with 
present value of the equation of state parameter w in the phantom regime. 
If this trend is verified by future datasets and if w remains in the phantom 
regime in the future then the fate of the universe acquires novel interesting 
features. The energy density of phantom fields increases with time and so does 
the predicted expansion acceleration rate ^. This monotonically increasing 
acceleration rate of the expansion may be shown to lead to a novel kind 
of singularity which occurs at a finite future time and is characterized by 
divergences of the scale factor a, the Hubble parameter H its derivative H 
and the scalar curvature. This singularity has been called 'Big Smash' jSHl 
the first time it was discussed and 'Big Rip' in a more recent study. 
An immediate consequence of the very rapid expansion rate as the Big Rip 
singularity is approached is the dissociation of bound systems due to the 
buildup of repulsive negative pressure in the interior of these systems. 
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This dissociation of bound systems can be studied by considering the 
spacetime in the vicinity of a point mass M placed in an expanding back- 
ground in order to study the effects of the cosmic expansion on bound sys- 
tems. Such a metric should interpolate between a static Schwarzschild metric 
at small distances from M and a time dependent Friedmann spacetime at 
large distances. In the Newtonian limit (weak field, low velocities) such an 
interpolating metric takes the form|52| : 

9G M 

ds 2 = (1 - . dt 2 - a{tf ■ {dp 2 + p 2 ■ {d6 2 + sin 2 0dtp 2 )) (45) 

a(t)p 

where p is the comoving radial coordinate. Using 

r = a(t) ■ p (46) 
the geodesies corresponding to the line element 11451 take the form 

a . GM 

_ (f _ _ r +rip 2 =0 (47) 

a r z 

and 

r 2 tp = L (48) 

where L is the constant angular momentum per unit mass. Therefore the 
radial equation of motion for a test particle in the Newtonian limit considered 
is 

d L 2 GM 



r = -r + — — (49) 

a r A 

The first term on the rhs proportional to the cosmic acceleration is a time 
dependent repulsive term which is increasing with time for w < — 1. This 
is easy to see by considering the Friedman equation H12|) combined with the 
dark energy evolution px ~ a~ 3 ( 1+tu ) where the scale factor obtained from 
the Friedman equation is 

a(t) = ^ — for t > t m (50) 

[-W + (1 + W)t/t m ] 3(1+™) 

and t m is the transition time from decelerating to accelerating expansion. For 
phantom energy (w < —1) the scale factor diverges at a finite time 

U = T^—t m > (51) 
1 + w 

leading to the Big Rip singularity. Clearly, the time dependent repulsive term 
of eq. iffify diverges at the Big Rip singularity. 

A quantitative analysis shows that the geodesic equation (|49|l is equiv- 
alent to a Newtonian equation with a time-dependent effective potential that 
determines the dynamics of the bound system which in dimesionless form is 

m 
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Veff = -^ + ^- ±A(t) V (52) 

where 

X(t) = ^±M_ (53 ) 

y ' 3(-w + (1 + w)t) v ; 



with w < — 1 and luq is defined as 



GM 



r o 

At t = 1 the system is assumed to be in circular orbit with radius given by 
the minimum r m i n (t) of the effective potential of equation l|52l) . It is easy to 
show that the minimum of the effective potential (|52|) disappears at a time 
t r i p which obeys 

16\/3Ty/2\l + 3w\ 
"9 6tt|1 



t* trip n n I I (55) 



The value of the bound system dissociation time t r i P may be verified by nu- 
merically solving the geodesic Newtonian equation of a test particle with the 
effective potential l|52[) . The resulting evolution close to the predicted dissoci- 
ation time tri p is shown in Fig. 16 for w — —1.2 and verifies the dissociation 
time predicted by equation (|55J) . Using the appropriate values for the bound 



Table 2. The difference between dissociation times t r i P and the big rip time t, for 
three bound systems in years as predicted by equation 1551 . The dissociation times 
t rip for the three bound systems in units of t m are also shown in column 3. The 
value w — —1.2 was assumed |53|. 



System 


t* — trip (yrs) 


trip 1 tm 


Solar System 


1.88 ■ 10 4 


6.00 


Milky Way 


3.59 ■ 10 s 


5.94 


Coma Cluster 


1.58 ■ 10 1U 


3.19 



system masses M the dissociation times of cosmological bound systems may 
be obtained. These are shown in Table 2. 



7 Future Prospects-Conclusion 

The question of the physical origin and dynamical evolution properties of dark 
energy is the central question currently in cosmology. Since the most sensitive 
and direct probes towards the answer of this question are distance-redshift 
surveys of Snla there has been intense activity during the recent years towards 
designing and implementing such projects using ground based and satellite 
observatories. Large arrays of CCDs such as MOSAIC camera at Cerro Tololo 
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Fig. 16. The numerically obtained evolution of a galactic size two body system at 
times close to the predicted dissociation time t r i P |58j . 

Inter- American Obsrevatory, the SUPRIME camera at Subaru or the MEGA- 
CAM at the Canada-France-Hawaii Telescope (CFHT) are some of the best 
ground based tools for supernova searches. These devices work well in the 
reddest bands (800-900nm) where the ultraviolet and visible light of redshifted 
high-z Snla is detected. Searches from the ground have the advantages of large 
telescope apertures (Subaru for example has 10 times the collecting area of 
the Hubble Space Telescope (HST)) and large CCD arrays (the CFHT has a 
378-milion pixel camera compared to the Advanced Camera for Surveys on 
HST which has 16 million pixels). On the other hand the advantage of space 
satellite observatories like the HST include avoiding the bright and variable 
night-sky encountered in the near infrared, the potential for much sharper 
imaging for point sources like supernovae to distinguish them from galaxies 
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in which they reside and better control over the observing conditions which 
need not factor in weather and moonlight. 

The original two Snla search teams (the Supernova Cosmology Project 
and the High-z Supernova Search Team) have evolved to a number of ongoing 
and proposed search projects both satellite and ground based. These projects 



\ Snla Projects 




SN Factory 
Carnegie SN Project 



ESSENCE 

CFHT Legacy Survey 



Higher-z SN Search 
(GOODS) 



SNAP 



Fig. 17. Ongoing and proposed Snla search projects with the corresponding redshift 
ranges. 



(see Fig. 17) include the following: 

• The Higher-z Supernova Search Team (HZT) [54] and the GOODS [55 
team of HST: This has originated from the High-z Supernova Search 
Team and has A. Riess of Space Telescope Sci. Inst, as its team leader. 
This team is in collaboration with the GOODS program (Great Obser- 
vatories Origin Deep Survey) using the ACS of the HST to detect and 
analyze high redshift (0.5 < z < 2) Snla. Successive GOODS observa- 
tions are spaced by 45 days providing 5 epochs of data on two fields: the 
Hubble Deep Field (HDF) north and south. Whereas the GOODS team 
adds these images to build a superdeep field, the HZT subtracts the ac- 
cumulated template image from each incoming frame. Thus the HZT has 
already detected more than 42 supernovae in the above redshift range. 

• Equation of State: SupEr Novae trace Cosmic Expansion (ESSENCE) [56 : 
This has also originated from the High-z Supernova Search Team and has 

C. Stubbs of the Univ. of Washington, C. Smith and N. Suntzeff of Cerro 
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Tololo as its team leaders. This ongoing program aims to find and measure 
200 Snla's in the redshift range of 0.15 < z < 0.7 where the transition from 
decelerating to accelerating expansion occurs. Spectroscopic backup to the 
program comes from the ground based Gemini, Magellan, VLT, Keck and 
MMT Obsevatory. The ESSENCE project is a five-year endeavor, with the 
goal of tightly constraining the time average of the equation-of-state pa- 
rameter w = p/ p of the dark energy. To help minimize systematic errors, 
all of their ground-based photometry is obtained with the same telescope 
and instrument. In 2003 the highest-redshift subset of ESSENCE super- 
novae was selected for detailed study with HST. 

The Supenova Legacy Survey (SNLS)HZ|: The CFHT Legacy Sur- 
vey aims at detecting and monitoring about 1000 supernovae in the red- 
shift range < z < 1 with Megaprime at the Canada-France-Hawaii 
telescope between 2003 and 2008. High-z spectroscopy of Snla is being 
carried on 8m class telescopes (Gemini, VLT, Keck). Team representa- 
tives are: C. Pritchet (Univ. Victoria), P. Astier (CNRS/IN2P3), S. Basa 
(CNRS/INSU) et. al. The SNLS has recently released the first year dataset 

M 

Nearby Supernova Factory (SNF)[5§|: The Nearby Supernova Fac- 
tory (SNF) is an international collaboration based at Lawrence Berkeley 
National Laboratory. Greg Aldering of Berkeley Lab's Physics Division is 
the principal investigator of the SNF. The goal of the SNF is to discover 
and carefully study 300 to 600 nearby Type la supernovae in the redshift 
range < z < 0.3. 

Carnegie SN Project (CSP)[59j: The goal of the project is the 
comprehensive study of both Type la and II Supernovae in the local 
(z < 0.07) universe. This is a long-term program with the goal of ob- 
taining exceedingly- well calibrated optical/near-infrared light curves and 
optical spectroscopy of over 200 Type la and Type II supernovae. The CSP 
takes advantage of the unique resources available at the Las Campanas 
Observatory (LCO). The team leader is R. Carlberg (Univ. of Toronto). 
Supernova Acceleration Probe (SNAP) [60 : This is a proposed space 
mission originating from LBNL's Supernova Cosmology Project that would 
increase the discovery rate for Snla's to about 2000 per year. The satellite 
called SNAP (Supernova / Acceleration Probe) would be a space based 
telescope with a one square degree field of view with 1 billion pixels. The 
project schedule would take approximately four years to construct and 
launch SNAP, and another three years of mission observations. SNAP has 
a 2 meter telescope with a large field of view: 600 times the sky area of the 
Hubble Space Telescopes Wide Field Camera. By repeatedly imaging 15 
square degrees of the sky, SNAP will accurately measure the energy spectra 
and brightness over time for over 2,000 Type la supernovae, discovering 
them just after they explode. 
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These projects aim at addressing important questions related to the physical 
origin and dynamical properties of dark energy. In particular these questions 
can be structured as follows: 

• Can the accelerating expansion be attributed to a dark energy ideal fluid 
with negative pressure or is it necessary to implement extensions of GR 
to understand the origin of the accelerating expansion? 

• Is w evolving with redshift and crossing the PDL? If the crossing of the 
PDL by w(z) is confirmed then it is quite likely that extensions of GR will 
be required to explain observations. 

• Is the cosmological constant consistent with data? If it remains consistent 
with future more detailed data then the theoretical efforts should be fo- 
cused on resolving the coincidence and the cosmological constant problems 
which may require anthropic principle arguments. 

The main points of this brief review may be summarized as follows: 

• Dark energy with negative pressure can explain Snla cosmological data 
indicating accelerating expansion of the universe. 

• The existence of a cosmological constant is consistent with Snla data but 
other evolving forms of dark energy crossing the w = — 1 line may provide 
better fits to some of the recent data (Gold dataset). 

• New observational projects are underway and are expected to lead to sig- 
nificant progress in the understanding of the properties of dark energy. 
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